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Since its introduction in the 1980s, magnetic resonance 
imaging (MRI) has become recognized as a powerful in 
vivo diagnostic tool. The objective of this article is to dis- 
cuss developments in quantitative MRI— and particularly 
texture analysis— that maximize diagnostic information. 
A fundamental part of the work involves careful study of 
the optimal MRI data collection strategies for texture 
analysis. This is critical, because different centers may vary 
their measuring sequences and acquisition protocols for 
clinical reasons, and may be reluctant to vary these for tex- 
ture investigation. Different measuring techniques, such 
as spin echo, gradient echo, and echo planar, and differ- 
ent measuring parameters produce totally different pat- 
terns in texture. Careful investigation of the dependence 
of all these variables using texture phantoms (test objects) 
will help understand how MRI image texture is formed 
from tissue structures. Therefore, it is essential to design 
and test reliable and accurate test objects for a detailed 
assessment of texture analysis methods in MRI. The main 
feature of these test objects is their ability to simulate tis- 
sue-like textures with tissue-like MR relaxation properties. 
Long-term stability is also vital, as is uniformity of the 
overall texture. Another aspect is to examine the test 
objects under a whole range of MRI measuring sequences 
and imaging conditions using different scanners to deter- 
mine their stability and utility. 
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n agnetic resonance imaging (MRI) is one of the 
most exciting imaging technologies for texture analysis: it 
offers the best soft tissue contrast, which can be dramati- 
cally varied during imaging. Careful study of the depen- 
dence of texture parameters on MRI data collection strat- 
egy is essential for texture analysis in order to avoid 
artificial texture from the scanner. This is critical, since dif- 
ferent centers may vary their measuring sequences and 
acquisition protocols for their clinical investigations. 
The basic problem in quantitative MRI texture analysis is 
the large number of different measuring techniques and 
imaging parameters, which can be easily changed during 
a clinical examination. Thus, different techniques and 
imaging parameters produce totally different patterns in 
the texture parameters of the same tissues in clinical 
examinations with different sensitivity to artificial texture 
overlaid by the scanner. The main problem in texture 
analysis with MRI is to avoid this artificial texture and 
minimize its influence. The presented work was per- 
formed in the framework of a European research project 
COST (Cooperation in the Field of Scientific and 
Technical Research) Bll between 1998 and 2002 by insti- 
tutions from 13 European countries, aimed at the devel- 
opment of quantitative methods for MRI texture analy- 
sis. 1 For further detail of texture analysis, parameters, and 
software, see the article by Materka in this volume 2 or ref- 
erences 3 to 7. 
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Selected abbreviations and acronyms 



EPI 


echo planar imaging 


FLASH 


fast low angle shot 


FOV 


field of view 


MA 


matrix size 


RF 


radio frequency 


ROI 


region of interest 


SNR 


signal-to-noise ratio 


TE 


spin echo time 


TH 


slice thickness 


TR 


repetition time 



Material and methods 

The complexity of this problem can be demonstrated by 
considering a typical measuring spin echo sequence as 
measured by a commercial whole-body imager. Various 
parameters can be easily changed during clinical investi- 
gation: image contrast is mainly defined by repetition time 
(TR) and spin echo time (TE); image resolution is defined 
by slice thickness (TH), field of view (FOV), and matrix 
size (MA), which also influence texture analysis. The para- 
meters of k-space acquisition and reconstruction are very 
important: k-space is the artificial space in which the raw 
MRI data are collected, and the image contrast and tex- 
ture is very sensitive to k-space strategies. Other parame- 
ters like coil setup and number of active coil segments are 
also responsible for signal and flip angle (a) variations in 
the image. Careful investigation of the dependence of all 
these variables will help understand how MRI image tex- 
ture is formed in tissue structures. In our studies, MRI 
acquisition was performed in the standard head coil of a 
1.5-T scanner (Siemens Vision, Erlangen, Germany). 

Spin echo technique 

One of the most important measuring techniques in clin- 
ical diagnosis is the spin echo sequence, in which 90° and 
180° radio frequency (RF) pulses produce the spin echo 
signal. In addition, gradients are used in x, y , and z direc- 
tions to localize the signal. 8 The advantages of this tech- 
nique are reduced artifacts, clearly defined contrast, and 
common availability. The disadvantages are the contrast 
dependency on RF pulse quality, and slice cross-talking, 
which is typical of a two-dimensional (2D) technique. 
This imaging technique allows measurement of the three 
relevant MRI tissue parameters: spin density (p), spin-lat- 
tice relaxation time (T^, and spin-spin relaxation time 



(T 2 ), which are most responsible for tissue contrast and 
texture. According to the theoretical equation for the spin 
echo signal 9 : 

S~p-(l-e TR/T1 )-e TEm [1] 

in which S is the spin echo signal, the contrast p can be cre- 
ated by a long TR and short TE, resulting in a flat image 
contrast and texture at high signal intensity (Figure la).Ti 
contrast can be created by short TR and short TE in spin 
echo imaging (Figure lb). On the other hand,T 2 contrast 
is created by long TR and long TE, mainly reflecting the 
water content of the tissue (Figure 1c). These three phys- 
ical tissue parameters are described in reference l.The real 
physical properties of tissues may be obscured by artificial 
contrast and texture from the scanner. 

Slice profile 

Slice profile is defined by the slice gradient and the 
shape of the RF pulse. Ideally, we would like to measure 
a rectangular slice, but due to technical reasons the real 
slice profile is Gaussian shaped. The consequence is that 
we have signal contributions from neighboring slices 
that influence the tissue texture. To minimize this effect, 
an interleaved slicing scheme is used in multislice 2D 
imaging. 

k-space 

Another aspect of artificial texture is connected to the k- 
space, which describes the strategy for raw data collec- 
tion. The k-space contains the measured signal frequen- 
cies k x and k y , the so-called hologram from which the real 
MRI image can be calculated by a Fourier transform 
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Figure 1. Spin echo images of a patient with meningioma. A. p-image 
(TR/TE = 2000 ms/1 0 ms). B. Ti image (TR/TE = 600 ms/1 0 ms). 
C. T 2 image (TR/TE 2000 ms/1 00 ms). TR, repetition time; TE, 
spin echo time. 
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(Figure 2). Some imaging techniques measure only every 
second line in the k-space to speed up the imaging 
sequence, which results in a reduction in the signal-to- 
noise ratio (SNR) by lHl and aliasing artifacts, with con- 
sequences for image texture. Restriction to the center of 
the k-space with zero filling of the outer part results in 
the same SNR effect without aliasing. 

RF excitation 

Another important variable is the RF characteristic and 
sensitivity of the transmitting and/or receiving coil, which 
can produce a lot of artificial texture from the scanner. 
This is demonstrated in Figure 3 using hard image scaling, 
which shows a clear signal inhomogeneity due to nonideal 
RF pulses at the outer range of the phantom (ie, coil). 
Another coil effect on image texture is produced with 
coil arrays, where the summarized image is a result of the 
combination of single coils, each of which contributes its 
own coil characteristics (eg, SNR, sensitivity, and RF exci- 
tation profile) to the summarized image. This means that 
image texture could slightly differ between the object 
center and the object boundary, where protons are close 
or far away from the center of the coil. 

Gradient echo techniques 

Significant effects on image contrast and texture are intro- 
duced by the imaging sequence itself, since the imaging sig- 
nal can have a very complex dependence on the physical 
properties of the underlying tissue. One example is the so- 
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Figure 2. The effect of k-space filling on image contrast and texture. The 
example demonstrates the strong dependence of image tex- 
ture on the k-space filling factor as used in so-called "keyhole" 
techniques. 



called gradient echo technique like FLASH (fast low angle 
shot), 10 where the 90° and 180° RF pulses are replaced by 
a low-angle RF pulse with a bipolar gradient scheme 
resulting in a gradient echo signal. This measuring tech- 
nique can be used as fast imaging 2D technique or as a real 
3D imaging technique because of the compact timing of 
the sequence. On the other hand, the FLASH signal has a 
complex dependence on T b the local spin-spin relaxation 
time (T 2 *), and the flip angle a, according to: 

S~A-e- TEm * 

I _ g-TR/Tl 

with A = p • since • [2] 

1- cosa - e~ TR/n 
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Figure 3. Effect of radio frequency (RF) profile of the head coil on image 
contrast and texture. A. Phantom measurement demonstrat- 
ing the signal inhomogeneity due to nonideal RF pulses at the 
outer range of the image (ie, coil). B to E. Result of a mea- 
surement using coil arrays, where the summarized image is 
a result of the combination of the single coils, which con- 
tribute with their own coil characteristic (signal-to-noise ratio, 
sensitivity, and RF excitation profile) to the contrast and tex- 
ture of the image. 



237 



Free paper 



Thus, the different flip angle distributions produced by 
the coil characteristics result in different signals, and as a 
consequence in different image texture patterns as 
demonstrated in Figure 4. 

A very complex signal and texture situation is present in 
so-called single shot imaging techniques like echo planar 
imaging (EPI), 11 where k-space is filled in one shot with 
multiple gradient echoes. This is achieved by a gradient 
scheme in which the upper corner of the k-space is 
reached by a single gradient pulse followed by a series of 
blips resulting in a rectangular movement through the k- 
space. This technique is very sensitive to local suscepti- 
bility artifacts, resulting in image distortions and strong 
T 2 * contrast dependence. 

Some special imaging techniques like spiral imaging can 
produce a very complex pattern in the image texture, 
since this single shot technique moves on a spiral through 
the k-space, which can be achieved by oscillating gradi- 
ents with a phase shift of 90° in the x and y directions. 
This technique requires data interpolation in k-space to 
bring the measured data onto a Cartesian coordinate sys- 
tem before Fourier transform. This interpolation can pro- 
duce spurious artifacts with the consequence that the 
image texture is dependent on k-space interpolation and 
image reconstruction. In addition, the problem of texture 
dependence on measuring technique is more complicated 



due to the large number of imaging sequences available 
on modern scanners, as illustrated in Figure 5. 

Results and discussion 

SNR dependence 

Figures 6a and 6b show the results of a FLASH experi- 
ment in a normal volunteer for SNR dependence mea- 
surement of texture parameters. The measuring parame- 
ters of the FLASH experiment were: TR/TE/a = 2 ms/ 
9 ms/30°; bandwidth (BW) = 195 Hz/pixel; MA = 
512x512; FOV = 280 mm;TH = 2 mm; and acquisitions 
(AQ) = 1 to 324 resulting in an SNR = 1 to 18. Texture 
parameters (SNR, entropy 5x5, correlation 5x5) of white 
matter, gray matter, and noise are shown as a function of 
the number of acquisitions (=SNR 2 ). Figure 6c demon- 
strates that no texture can be measured in white matter 
using standard image resolution (0.5x0.5x2 mm 3 ) as 
described above, since the SNR of white matter has the 
same characteristics as noise. In contrast, the SNR of gray 
matter reaches a nearly constant value at about 16 acqui- 
sitions and no further improvement can be reached due 
to the true underlying texture of the tissue. The same 
observation holds for a typical parameter of microtex- 
ture, like entropy 5x5 (Figure 6d), while no dependence 




Figure 4. Example of a 3D FLASH (fast low angle shot) dataset of a normal volunteer measured by the conventional head coil. A. Excellent Ti con- 
trast in the original transversal images. B. A strong signal inhomogeneity is obvious in the reconstructed sagittal images as seen in the 
homogeneous phantom and corresponding anatomical structures, like upper part of the brain and cerebellum. This artifact is due to the 
radio frequency excitation profile of the head coil, which produces different Ti weightings and signal amplitudes as a function of the 
flip angle according to the FLASH formula (Equation 2). 
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on SNR can be detected for a typical parameter of 
macrotexture, like correlation 5x5 (Figure 6e). Based on 
this observation, a sufficient SNR>4 is necessary to mea- 
sure the real textural behavior of the human brain. 12 13 

Normalization 

A texture test object (PS AG) was developed on the basis 
of polystyrene (PS) and agar solution (AG) to mimic tex- 
ture properties artificially. PS spheres are available from 
the technological process of PS production. Two types of 
spheres were used for the phantom construction: ran- 
domly distributed spheres of diameter 0.2 to 3.15 mm; or 
mechanically separated spheres of diameter 0.8 to 1.25 
mm, 1.25 to 2 mm, or 2 to 3.15 mm. Polyethylene tubes of 
diameter 1.5 and 2.8 cm were filled with spheres and by 




a hot solution of 4% agar (free and doped with DyCl 3 ). 
One milliliter of 0.1% NaN 3 was added per liter of agar 
for microbiological stability. 14 

A second texture test object containing foam at different 
densities in Gd-DTPA solution was used to describe 
microtexture properties. Phantom tubes containing foams 
with coarse, middle, and fine density were constructed and 
filled with a Magnevist® (Schering, Berlin) solution at a 
concentration of 1:4000. Problems with the foam phantoms 
are air bubbles, which create susceptibility artifacts in the 
images, and so a careful preparation of the foam phantoms 
is necessary. Both types of phantoms were placed next to 
the head of a volunteer and a position for the imaging slab 
was chosen such that all vials and part of the volunteer's 
brain were contained in the 3D slab. With this setup sev- 
eral 3D data sets with different imaging parameters were 




Figure 5. Sketch of the family of imaging techniques available on modern scanners. There are many strategies of mixing spin echoes with gradient 
echoes to speed up imaging time with the consequence of very complex image contrast and texture. CISS, constructive interference in the 
steady state; CSE, conventional SE; DESS, dual echo steady state; EPI, echo planar imaging; FAST, Fourier-acquired steady state; FISP, fast imag- 
ing with steady precession; FLASH, fast low angle shot; FSE, fast spin echo; FSPGR, fast spoiled gradient recalled acquisition into steady state; 
GRASE, gradient and spin echo; GRASS, gradient recalled acquisition into steady state; GSE, gradient SE; HASTE, half-Fourier acquisition 
single-shot turbo SE; IR, inversion recovery; MPRAGE, magnetization prepared gradient echo imaging; RARE, rapid acquisition with relaxation 
enhancement; SE, spin echo; SPGR, spoiled gradient recalled acquisition into steady state; Ti-FFE, Ti-weighted fast field echo; T2-FFE, T2- 
weighted fast field echo; TFLAIR, turbo fluid attenuated IR; TGSE, turbo GSE; TIR, turbo IR; TIRM, turbo IR magnitidue; TSE, turbo SE. 
Copyright © Siemens. 
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acquired to demonstrate the influence of resolution and 
SNR, as well as the dependence of the texture parameters 
on different imaging parameters (eg, oc,TR,TE). In a pilot 
study, texture parameters such as mean gradient show the 
same behavior in phantoms as in white matter for differ- 
ent patients, indicating that a normalization of texture 
parameters using test objects is possible (Figure 7). 
However, texture normalization is necessary, but it is not 
possible to mimic all texture features by phantoms. 15 
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Figure 6. FLASH (fast low angle shot) images of a normal volunteer for 
measuring signal-to-noise (SNR) dependence of texture para- 
meters at (A) SNR =1 (1 acquisition) and (B) SNR =1 8 (324 acqui- 
sitions). C to E. Texture parameters (SNR, entropy 5x5, correla- 
tion 5x5) of white matter, gray matter, and noise are shown as 
a function of the number of acquisitions (=SNR 2 ). SD, standard 
deviation. 



Figure 7. A. Three-dimensional FLASH (fast low angle shot) image of a 
patient with glioblastoma with texture test objects located 
beside the head for testing texture normalization. B. Texture 
parameters such as mean gradient show the same behavior in 
the phantom as in white matter for different patients, indi- 
cating that normalization of texture parameters using test 
objects is possible. 
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Clinical application 

The aim of this pilot study was to assess the possibility of 
quantitative description of texture directivity in trabecu- 
lar bone with an attempt to quantitative description of 
trabecular bone structural anisotropy using texture 
analysis of 3D FLASH MRI. A series of 3D FLASH 
images, all of 256x256 pixels, with the voxel size of 
0.4x0.4x0.4 mm 3 , were measured on a standard 1.5-T 
scanner (Siemens Vision, Erlangen, Germany) using a 
small flex coil. The images in Figure 8 represent trabecu- 
lar bone cross-sections in the sagittal and reconstructed 
transversal direction. For bone image texture analysis, cir- 
cular regions of interest (ROI) were marked on corre- 
sponding bone cross-sections and effort has been made 




Horizontal Vertical 45° 135° 



Figure 8. Example of a trabecular bone 3D FLASH (fast low angle shot) 
measurement of a normal volunteer in sagittal (A) and recon- 
structed transversal direction (B) for testing texture directivity. 
(C) Texture parameters like run length nonuniformity show a 
clear dependence on direction of gravity. Quantitative analy- 
sis of this directivity is important to medical diagnosis, eg, in 
early detection of osteoporosis, as the directivity may vary 
according to the development of the disease. 



to maintain a large-size ROI for better statistical signif- 
icance of texture parameters. The texture of the bone 
image shows apparent directivity, which reflects 
anisotropy of its physical structure according to the direc- 
tion of gravity (Figure 8c). Quantitative analysis of this 
directivity is important to medical diagnosis, eg, in early 
detection of osteoporosis, as the directivity may vary 
according to the development of the disease. □ 
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Problemas en el analisis de la textura con ima- 
genes de resonancia magnetica 



Las imagenes de resonancia magnetica (IRM) se han 
reconocido como una poderosa herramienta para 
el diagnostico in vivo desde su introduccion en la 
decada de 1980. El objetivo de este artfculo es dis- 
cutir el desarrollo de las IRM de tipo cuantitativo y 
en particular el analisis de la textura, lo que permite 
maximizar la informacion para el diagnostico. Una 
parte fundamental de este trabajo incluye un cui- 
dadoso estudio de las mejores estrategias de reco- 
leccion de datos de las IRM para el analisis de la tex- 
tura. Este aspecto es central ya que diferentes 
centros, por razones clmicas, pueden variar la 
secuencia de medicion y los protocolos de adquisi- 
cion de datos y por lo tanto pueden ser reacios a 
modificar estos para la investigacion de la textura. 
Diversas tecnicas de medicion como el eco de spin, 
el eco de gradiente y el eco planar ydistintos para- 
metros de medicion producen patrones de textura 
totalmente diferentes. Una investigacion cuidadosa 
de la influencia de todas estas variables mediante 
el empleo de texturas fantasma (objetos de prueba) 
ayudara a comprender como se forman las texturas 
de las IRM a partir de las estructuras de los tejidos. 
Por lo tanto f resulta esencial disehar y ensayar obje- 
tos de prueba exactos y confiables que permitan 
una evaluacion detallada de los metodos de anali- 
sis de la textura de las IRM. La caractenstica princi- 
pal de estos objetos de prueba es su capacidad para 
simular texturas parecidas a los tejidos y con pro- 
piedades semejantes a estos ultimos. La estabilidad 
a largo plazo de estos modelos es de gran impor- 
tancia, como tambien la uniformidad global de la 
textura. Otro aspecto consiste en examinar la esta- 
bilidad y utilidad de los objetos de prueba 
mediante diversas secuencias de medicion y condi- 
ciones en que se realicen las IRM empleando dis- 
tintos resonadores. 



Problemes poses par les analyses de texture 
au cours de I'imagerie par resonance 
magnetique 

Depuis son introduction au milieu des annees 
quatre-vingt, I'imagerie par resonance magnetique 
(IRM) est devenue un puissant outil de diagnostic in 
vivo. L'objectif de cet article est d'evaluer les avan- 
cees de I'IRM quantitative, notamment en ce qui 
concerne les analyses de texture, qui optimisent ^in- 
formation diagnostique. Cet article sera en grande 
partie consacre a I'etude minutieuse des meilleures 
strategies de recueil de donnees d'IRM pour les ana- 
lyses de texture. Ceci est essentiel, dans la mesure ou 
les sequences de mesures et les protocoles de saisie 
cliniques peuvent varierselon les centres et que ces 
derniers peuvent etre reticents a les modifier pour 
les analyses de texture. Des differences aussi bien 
dans les techniques de mesure utilisees, qu'il s'agisse 
de I'echo de spin, I'echo de gradient, I'echo planaire, 
etc., que dans les parametres mesures retentissent 
sur la texture en entrainant /'apparition de motifs 
totalement differents. devaluation precise de I'in- 
fluence de ces variables en utilisant des images eta- 
Ion (« fantdmes ») de texture (ou objets de test) faci- 
litera une meilleure comprehension du mode de 
constitution des textures des structures tissulaires en 
IRM. Par consequent, il est essentiel de developper 
et d'experimenter des objets de test exacts et 
fiables pour une evaluation detaillee des methodes 
utilisees en IRM pour les analyses de texture. La prin- 
cipale caracteristique de ces modeles reside dans leur 
capacite a imiter des textures semblables a celles pro- 
duces par les tissus et ayant des proprietes de relaxa- 
tion RM identiques. La stabilite a long terme de ces 
modeles revet egalement une grande importance, 
tout comme I'uniformite globale des textures qui en 
resultent. Un autre objectif est d'examiner la stabi- 
lite et I'utilite de ces objets de test en les soumettant 
a toute une gamme de sequences de mesures et de 
conditions d'imagerie par IRM et de les etudier avec 
divers types d'appareils d'IRM. 
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